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Abstract. The use of femtosecond laser pulses generated by a Ti: Sapphire laser 
system allows us to gain an insight into the magnetization dynamics on time scales from 
sub-picosecond up to 1 ns directly in the time domain. This experimental technique is 
used to excite a polycrystalline nickel (Ni) film optically and probe the dynamics 
afterwards. Different spin wave modes (the Kittel mode, perpendicular standing 
spin- wave modes (PSSW) and dipolar spin- wave modes (Damon-Eshbach modes)) are 
identified as the Ni thickness is increased. The Kittel mode allows determination of the 
Gilbert damping parameter a extracted from the magnetization relaxation time Tq,. 
The non-local damping by spin currents emitted into a non-magnetic metallic layer 
of vanadium (V), palladium (Pd) and the rare earth dysprosium (Dy) are studied 
for wedge-shaped Ni films 1 nm — 30 nm. The damping parameter increases from 
a = 0.045 intrinsic for nickel to a > 0.10 for the heavy materials, such as Pd and 
Dy, for the thinnest Ni films below 10 nm thickness. Also, for the thinnest reference Ni 
film thickness, an increased magnetic damping below 4nm is observed. The origin of 
this increase is discussed within the framework of line broadening by locally different 
precessional frequencies within the laser spot region. 
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1. Introduction 

The understanding of picosecond-pulsed excitation of spin packets, spin wave modes 
and spin currents is of importance in developing a controlled magnetic switching concept 
beyond the hundred picosecond timescale and to test the speed of magnetic data storage 
media heading to the physical limits. Over the last years profound progress has been 
made within that field by using femtosecond laser spectroscopy. The recent discoveries 
in ultrafast magnetization dynamics are heading to a new understanding [15] and 
new all-optical switching concepts have been discovered [6]. In addition, the all-optical 
method has developed into a valuable tool to study the magnetization dynamics of 
the magnetic precession and thereby access magnetocrystalline anisotropics and the 
magnetic damping [7-11] or the dynamics of magnetic modes in nanometer sized arrays 
of magnetic structures [12,13] and single magnetic nanostructures [14,15]. Naturally, 
one finds similarities and differences as compared to magnetic resonance techniques 
in frequency space (FMR) [16], optical techniques such as Brillouin light scattering 
(BLS) [17,18] and time-resolved techniques, for example pulsed inductive magnetometry 
(PIMM) [19]. Advantages and disadvantages of the different techniques have already 
been compared in previous work [20 22]. The same concepts can be applied to the 
femtosecond-laser-based all-optical spectroscopy techniques. Here we discuss their 
abilities, highlighting some aspects and peculiarities [11,23-27]: 

i. After excitation within the intense laser pulse, the nature of the magnetic relaxation 
mechanisms determine the magnetic modes observed on the larger time scale [5]. 
For a Ni wedge different modes are found as the thickness is increased: coherent 
precession (Kittel mode), standing spin waves (already found in [28]) and dipolar 
surface spin waves (Damon-Eshbach modes) appear and can be identified. 

ii. Magnetic damping has been extracted by the use of fs spectroscopy experiments 
already in various materials, epitaxial films, as a function of the applied field 
strength, field orientation and laser excitation power [7-11]. Using the Kittel 
mode, we study the energy dissipation process caused by non-local damping by spin 
currents [29] in Ni by attaching a transition metal film (vanadium (V), palladium 
(Pd) and a rare earth film (dysprosium (Dy)) as a spin sink material and compare 
them to a Ni reference sample. The present advantages and disadvantages of the 
method are discussed. 

iii. A modification of the magnetic damping is found for the thinnest magnetic layers 
below 4nm. The understanding of this effect is of high interest because of the 
increase in methods used to study magnetic damping processes in the low field 
region in the current literature. We present a simple model of line broadening 
known from FMR [30-32] and adapted to the all-optical geometry that pictures 
the effect of the increased intrinsic apparent damping observed. Therein a spread 
local magnetic property within the probe spot region is used to mimic the increased 
apparent damping for the low field region. 
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Figure 1. a) Schematics of the pump probe experiment to determine the change in 
Kerr rotation as a function of the delay time r. b) Experimental data on short and 
long time scales. On top a schematic on the processes involved is given. 



2. Experimental Technique 

The all-optical approach to measuring magnetization dynamics uses femtosecond laser 
pulses in a pump-probe geometry. In our experimental setup a Ti: Sapphire oscillator 
generates the fs laser pulses which are then amplified by a regenerative amplifier (RegA 
9050). This systems laser pulse characteristics are 815 nm central wave length, a 
repetition rate of 250 kHz, a temporal length of 50 — 80 fs and an energy of ^ 1 /iJ 
per pulse. The beam is split into a strong pump beam (95% of the incoming power), 
which triggers the magnetization dynamics by depositing energy within the spot region. 
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and a weaker probe pulse (5% of the incoming power) to probe the magnetization 
dynamics via the magneto-optical Kerr effect delayed by the time r, in the following 
abbreviated as time-resolved magneto-optic Kerr effect (TRMOKE). The schematic 
setup and sample geometry is given in figure [T^) . The spot diameters of the pump 
and probe beam are 60 /xm and 30 /xm respectively. A double-modulation technique is 
applied to detect the measured signal adapted from [33] : the probe beam is modulated 
with a photo-elastic modulator (PEM) at a frequency /i = 250 kHz and the pump 
beam by a mechanic chopper at a frequency /2 = 800 Hz. The sample is situated in 
a variable magnetic field (0 — 150 mT), which can be rotated from 0° (in-plane) to 90° 
(out of plane) direction. The degree of demagnetization can be varied by the pump 
fluence (lOmJ/cm^ — 60mJ/cm^) to up to 20% for layer thicknesses around 30 nm and 
up to over 80% for layers thinner than 5nm. The samples studied were all grown on 
Si(lOO) substrates by e-beam evaporation in a UHV chamber at a base pressure of 
^ 5 X 10~^° mbar. For a variation of the thickness, the layers are grown as wedges with 
a constant gradient on a total wedge length of 15 mm. 

3. Results and discussion 

3.1. Kittel mode^ standing spin waves and Damon- Eshbach surface modes 

To give an introduction to the TRMOKE signals AOxerri^) measured on the timescale 
from picoseconds to nanoseconds first, the ultrafast demagnetization on a characteristic 
time scale tm and the magnetic precessional motion damped on a time scale is shown 
for a Ni film in figure ^p) ; the schematics of the processes involved on the different time 
scales are given on the top. The change in Kerr rotation AOxerri^) shows a sudden drop 
at r = ps. This mirrors the demagnetization within a timescale of ^ 200 fs [34-36] . For 
the short time scale the dynamics are dominated by electronic relaxation processes, as 
described phenomenologically in the three temperature model [34] or by connecting the 
electron-spin scattering channel with Elliot- Yafet processes, as done by Koopmans [36] 
and Chantrell [4] later. At that time scale the collective precessional motion lasting up 
to the nanosecond scale is initiated [28,37]: the energy deposited by the pump pulse 
leads to a change in the magnetic anisotropy and magnetization, and thus the total 
effective field. Within ^ 10 ps the total effective field has recovered to the old value 
and direction again. However, the magnetization, which followed the effective field, 
is still out of equilibrium and starts to relax by precessing around the effective field. 
This mechanism can be imagined as a magnetic field pulse a few picoseconds long, and 
is therefore sometimes called an anisotropy field pulse. The resulting anisotropy field 
pulse is significantly shorter than standard field pulses [38] . This makes the TRMOKE 
experiment different to other magnetization dynamics experiments. 

The fact that the situation is not fully described by the model can be seen in the 
following. Already van Kampen et al. [28] not only observed the coherent precessional 
mode, they also identified another mode at a higher frequency than the coherent 
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precession mode, shifted by ujk,n ^ 2Ak'^ = 2Ami /t^-^ ^ the standing spin wave (PSSW) 
mode. It originates from the confinement of the finite layer thickness, where A is the 
exchange couphng constant and n is a given order. Here we also present the finding of 
dipolar propagating spin waves. For all three, the frequency dependence as a function 
of the applied magnetic field will be discussed, a necessity for identifying them in the 
experiments later on. 

For the coherent precession the frequency dependence is described by the Kittel 
equation. It is derived by expressing the eflFective field in the Landau-Lifshitz-Gilbert 
(LLG equation) as a partial derivative of the free magnetic energy [39,40]. Assuming 
negligible in-plane anisotropy in case of the polycrystalline nickel (Ni) film and small 
tilting angles of the magnetization out of the sample plane (field is applied 35° out of 
plane figure [T^)), it is solved as derived in [41]: 



^ ^ Mo V ^^'^ (^Moi/x + /XqM, - (1) 

For the standing spin waves (PSSW) a similar equation is given. For the geometry 
with the field applied 35° out of plane (figure [T^) the frequencies uj and ujk.n do not 
simply add as in the field applied in plane geometry [41]: 



7 /, 2Ak^^f _ 2K, 2Ak^\ 

. = l.^(,„H^ + —) + - m7 + m7 j . (2) 

While the exchange energy dominates in the limit of small length scales, the 
magnetic dipolar interaction becomes important at larger length scales. Damon and 
Eshbach [42] derived by taking into account the dipolar interactions in the limit of 
negligible exchange energy, the solution of the Damon-Eshbach (DE) surface waves 
propagating with a wave vector q along the surface, decaying within the magnetic layer. 
The wavelengths are found to be above the > jim range for Ni [27] . 




^ = — i/mo^x ( Mo^x + I^^Ms - -j^ + ^ [1 - exp{-2qtNi)]U (3) 

The depth of the demagnetization by the femtosecond laser pulse is given by the 
optical penetration length Aopt ^ 15 nm (A = 800 nm). From the nature of the excitation 
process in the TRMOKE experiment one can derive that for different thicknesses ^Ni it 
will change from an excitation of the full film for a ^ 10 nm film to a thin excitation 
layer only for a few 100 nm thick film; thus the excitation will be highly asymmetric. 
The model of the magnetic anisotropy field pulse fails to explain these eflFects since it is 
based on a macrospin picture. 

Another way to look at the excitation mechanism has been discussed by Djordjevic 
et al. [5]. When the magnetic system is excited, on a length scale of the optical 
penetration depth short wavelength (high k vector) spin- wave excitations appear. As 
time evolves, two processes appear: the modes with high frequency owning a fast 
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oscillation in space are damped very fast by giving part of the deposited energy to 
the lattice. In addition, through multiple magnon interaction lower fc-vector states are 
populated, resulting in the highest occupation of the lowest energy modes at the end 
(e.g. the PSSW and DE modes here). As the Ni thickness is increased, the excitation 
profile becomes increasingly asymmetric, favoring inhomogeneous magnetic excitations, 
as the PSSW mode. The DE modes, due to their nature based on a dipolar interaction, 
are expected to be found only for higher thicknesses. 
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Figure 2. Change in Kerr rotation after excitation on the long time scale for Cu 2nm/ 
NitNi nm/Si(100) with ^Ni = 20 and b) their Fourier transform for different applied 
fields — 150 mT, (35° out of plane (blue)). In c) the Fourier power spectra as color 
maps for three Ni thicknesses ^Ni = 20, 40 and 220 nm are given. The data overlaid is 
determined form the peak positions. The straight lines are the analysis of the different 
modes and are identified in the graph (Kittel model), perpendicular standing spin wave 
(PSSW) and dipolar surface spin wave (Damon Eshbach mode). 
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The identification of the mode is important in determining a value for the magnetic 
damping a. Figure [2] pictures the identification of the different modes and their 
appearance for diflFerent Ni thicknesses. The data are handled as follows: for a 
t^i = 20 nm film on Si(lOO), covered with a 2nm Cu protection layer, in a) the original 
data after background subtraction and in b) its corresponding Fourier transform, shown 
for increasing applied magnetic field. The evolution of the mode frequency and its 
amplitude increase can be followed. An exponentially decaying incoherent background 
is subtracted from the data. This has to be done very carefully, to avoid a step- 
like background which will be evident after Fourier transform as a sum of odd higher 
harmonics. The frequency resolution is limited by the scan range of 1 ns corresponding 
to A(jj/27r = 1 GHz. However, since the oscillation is damped within the scan range, 
the datasets have been extended before Fourier transform to increase their grid points. 
A color map of the power spectrum is shown in figure [2]3), where the peak positions are 
marked by the data points overlaid. For the 20 nm thick film with t^i = 20 nm ~ Aopt 
only a single mode is observed. The mode is analyzed by [T] indicating the Kittel mode 
being present (data points and line in figure [2]3), top) using = 3.03 • 10^ J/m^. With 
increasing nickel thickness t^i = 40 nm > Aopt, the perpendicular standing spin waves 
(PSSW) of first order are additionally excited and start to appear in the spectra (figure 

middle). An exchange constant A = 9.5 • 10^^ J/m is extracted. In the limit of 
t^i = 220 nm ^ Aopt (figure [l]^)) the excitation involves the surface only. Hence, modes 
with comparable amplitude profile, e.g. with their amplitude decaying into the Ni layer, 
are preferred. Consequently DE surface waves are identified as described by |3] and 
dominate the spectra up to critical fields as high as /xoi^crit = 100 mT. For t^i = 220 nm 
the wave factor is /c = 2 /xm (data points and line in figure [2]3), bottom). For larger 
fields than 100 mT the DE mode frequency branch merges into the Kittel mode [27]. 

To resume the previous findings for the first subsection, we have shown that in 
fact the DE modes, though they are propagating spin- wave modes, can be identified 
in the spectra and play a very important role for Ni thicknesses above ^Ni = 80 nm. 
They appear for thicknesses much thinner than the wavelength of the propagating 
mode. Perpendicular standing spin waves (PSSW) give an important contribution to 
the spectra for Ni thicknesses above t^i = 20 nm. For thicknesses below t^[ = 20 nm we 
observe the homogeneously precessing Kittel mode only. This thickness range should 
be used to determine the magnetic damping in TRMOKE experiments. 

3.2. Data analysis: determination of the magnetic damping 

For the experiments carried out in the following with t^i < 25 nm the observed dynamics 
can be ascribed to the coherent precession of the magnetization (Kittel mode). The 
analysis procedure is illustrated in the following using the data given in figure [3^). A 
Pd layer is attached to a Ni film with the thicknesses (NilOnm/ Pd5nm/ Si(lOO)) to 
study the non-local damping by spin currents absorbed by the Pd. The different spectra 
with varying the magnetic field strength from OmT — 150 mT are plotted from bottom 
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Figure 3. a) Kerr rotation spectra for a Cu2nm/ NilOnm/ Pd 5 nm/Si(100) layer, 
measured for fields applied from — 150 mT (35° out of plane, (blue)) and the fitted 
functions (white, dashed), b) The magnetic damping a and precession frequencies 
extracted from the fits to the measured spectra. The line is given by the Kittel mode 
(gray). 

The data can be analyzed using the harmonic function with an exponential decay 
within Ta'. 

A9k - exp (^-^^ ' sin(27r(r - To)iy) + 5(r), (4) 

The precession frequency u = uj /27r and the exponential decay time of the 
precession amplitude is extracted, where the function B{r) stands for the background 
arising from the uncorrelated magnetic and phonon excitations. To determine the 
Gilbert damping parameter a as given in the ansatz by Gilbert, the exponential decay 
time Ta has to be related with a. The LLG equation is solved under the same 
preconditions as for equation [T] using an exponential decay of the harmonic precession 
within Ta from|4} Then the damping parameter a and can be expressed by the following 
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equation [41]: 



It is evident from [5] that in order to determine the Gilbert damping a from the decay 
of the Kittel mode r^, the variables 7, Mg and have to be inserted, and therefore 
has to be determined beforehand. 

In figure [3^), the background B{t) is already subtracted. The fits using [i] are 
plotted with the dashed lines on top of the measured spectra. The results are presented 
in b). The frequencies range from 3 GHz for 30 mT to 7.5 GHz for the 150 mT applied 
magnetic field. They increase linearly with the strength of the applied magnetic field for 
high field values. The extrapolated intersection with the ordinate is related to the square 
root of the dipolar and anisotropy field. Using the Kittel equation ([T]), one determines 
the out-of-plane anisotropy constant of = 6.8- 10^ J /m^. The calculated magnetic 
damping a as a function of the applied field is given in the graph below: this is mostly 
constant but increases below 60 mT. Within the ansatz given by Gilbert, the damping 
constant a is assumed to be field-independent. We find that this is fulfilled for most 
of the values: the average value of a = 0.0453(4), consistent with earlier findings by 
Bhagat and Lubitz from FMR experiments [43], is indicated by the line in the plot. The 
a given in the following will always be averaged over a field region where the damping 
is Gilbert-like. A deviation from this value occurs for the small external field strengths. 
It originates for two reasons: fitting [5] with a few periods only does not determine a 
reliable value of the exponential precession decay time and leads to a larger error. 
Second, magnetic inhomogeneities mapping a spread in anisotropy energies within the 
probe spot region can also be a source, and this becomes generally more important for 
even thicker films below 4nm [31]. This will be discussed in more detail in the last 
section of the manuscript. 

3.3. Intrinsic damping: nickel wedge 

For our experiments Ni was chosen instead of Fe or Py as a ferromagnetic layer. The 
latter would be preferable because of their lower intrinsic damping aint, which make 
the films more sensitive for detecting the non-local contribution to the damping. The 
reason for using Ni for our experiments is the larger signal excited in the TRMOKE 
experiments. The magnetic damping a[^t is used as a reference later on. The different 
spectra with varying the Ni thickness t^[ Nixnm/Si(100) from 2nm < x < 22 nm are 
plotted from bottom to top (with the constant magnetic field 150 mT and tilted 30° 
out of plane) in figure |4^) . The measurements were performed immediately after the 
sample preparation, in order to prevent oxidation on the nickel surface caused by the 
lack of a protection layer (omitted on purpose). The spectra show similar precession 
frequency and initial excitation amplitude. However, the layers with t^i < 10 nm show a 
frequency shift visually recognized in the TRMOKE data. Furthermore, the precession 
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amplitude decreases faster for the thinner layers. Figure [5] shows the frequencies and 
the damping parameter extracted from the measured data in the intrinsic case for the 
nickel wedge sample (black squares). While the precession frequency given for 150 mT 
is almost constant above 8nm Ni thickness, it starts to drop by about 25% for the 
thinnest layer. The magnetic damping a (black squares) is found to increase to up to 
a = 0.1, an indication that in addition to the intrinsic there are also extrinsic processes 
contributing. It has to be noted that the change in a is not correlated with the decrease 
of the precession frequency. The magnetic damping a is found to increase below a 
thickness of 4nm, while the frequency decrease is observed below a thickness of 10 nm. 
A priori 7, Mg and can be involved in the observed frequency shift, but they can 
not be disentangled within a fit of our field-dependent experiments. However, from our 
magnetic characterization no evidence of a change of 7 and Mg is found. A saturation 
magnetization fi^Mg = 0.659 T and ^'-f actor of 2.21 for Ni are used throughout the 
manuscript[|] and is determined as a function of the Ni thickness, which shows a 
l/^Ni behavior, as expected for a magnetic interface anisotropy term [44]. 

The knowledge of the intrinsic Gilbert damping ai^t of the Ni film of a constant 
value for up to 3 nm thickness allows us to make a comparative study of the non-local 
damping a^, introduced by an adjacent layer of vanadium (V) and palladium (Pd) as 
representatives for transition metals, and dysprosium (Dy) as a representative of the 
rare earths. Both damping contributions due to intrinsic a[^t and non-local spin current 
damping are superimposed by: 

= OLint + ol'- (6) 

They have to be disentangled by a study of the thickness dependence and compared 
to the theory of spin-current pumping, plus a careful comparison to the intrinsic value 
Qint has to be made. 

3.4' Non-local spin current damping: theory 

Dynamic spin currents excited by a precessing moment in an adjacent nonmagnetic 
layer (NM) are the consequence of the fact that static spin polarization at the interface 
follows a dynamic movement of a collective magnetic excitation. The eflFect has already 
been proposed in the seventies [45,46] and later calculated within a spin reservoir model 
with the spins pumped through the interfaces of the material by Tserkovnyak [29,47]. 
For each precession, pumping of the spin current results in a corresponding loss in 
magnetization, and thus in a loss of angular momentum. The spin information is lost 
and the backward diflFusion damps the precession of the magnetic moment. In addition to 
the first experiments using ferromagnetic resonance (FMR) [48-54] it has been observed 
in time- resolved experiments using magnetic field pulses for excitation [55,56]. In fact 

X An altered g-factor by interface intermixing can not decrease its value below ~2. Also, there is no 
evidence for a reduced Ms for lower thicknesses found in the Kerr rotation versus Ni thickness data. 
More expected is a change in the magnetic anisotropy Kz. For the calculation of a later on, the in 
both cases (assuming a variation of Kz or an altered 7) the differences are negligible. 
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Figure 4. a) Kerr rotation spectra for nickel layers from ^Ni = 2 nm — 22 nm, measured 
on the nickel wedge ^Ni = nmNi/Si(100) and opposed in b) by a nickel wedge A12nm/ 
DySnm/ t^i = nmNi/Si(100) with a 5nm Dy spin-sink layer. 



the non-local spin current damping is very closely related to the damping by spin-flip 
scattering described within the s-d current model [57, 58] that uses the approximation 
of strongly localized d-states and delocalized s-states [59]. 

A review describes the underlying circuit theory and dynamics of the spin currents 
at interfaces in detail [60]. The outcome of the theoretical understanding is that the 
additional Gilbert damping is proportional to the angular momentum A^^i transmitted 
through the interface. Since each interface owns a characteristic reflection and 
transmission, the size of A^^i depends on the matching of the Fermi surfaces. The 
absolute value is given by the total balance between transmitted angular momentum 
and the back flow. For the non-local damping one finds: 



1 



a 



47rM,tFM 



1 



^^tanh 



(7) 



The tanh function stems thereby from the diffusion profile of the spin currents 
determined by the spin diffusion length Agd within the non-magnetic material with 
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thickness ^nm- Also, one finds from the analysis the ratio of the electron scattering 
rate Tgi versus the spin fiip rate Tgf. The total amount of spin current through the 
interfaces is determined by the interface spin mixing conductance G^K It is related 
to the magnetic volume. It is therefore that scales with the thickness of the magnetic 
layer tpM- The effective gyromagnetic ratio altered by the spin-current implies that in 
addition to an increased damping a small frequency shift will be observed. The non-local 
Gilbert damping becomes important when it exceeds the intrinsic damping aint- 

3.5. Non-local damping: vanadium, palladium and dysprosium 

Differing from other techniques, TRMOKE experiments require optical access for 
excitation and detection, setting some restrictions to the layer stack assembly that can 
be investigated with this method: a thick metallic layer on top of the magnetic layer is 
not practical. Placing the damping layer below the magnetic layer is also unfavorable: 
by increasing the spin sink thickness the roughness of the metal film will increase with 
the metals layer thickness and introduce a different defects density, altering aint- In the 
following the nickel thickness will be varied and the spin sink thickness will be kept fixed 
at 5nm. To warrant that the nickel films magnetic properties are always comparable 
to the reference experiment (X^, o^mt), they are always grown first on the Si(lOO). For 
the Pd case the damping layer is below the Ni layer. Here the excitation mechanism 
did not work and the oscillations were too weak in amplitude to analyze the damping 
a, probably due to the high refiectivity of Pd. 

The results are presented in figure|4]3) for the nickel wedge sample Nixnm/Si(100) 
with a 5 nm dysprosium (Dy) as a spin sink layer, covered by an aluminum protection 
layer, as opposed to the nickel wedge sample data without this in a). The nickel layer 
thickness is varied from 2 nm < x < 22 nm. All spectra were measured in an external 
magnetic field set to 150 mT and tilted 30° out of plane. For the thinnest Ni thickness, 
the amplitude of the precession is found to be smaller due to the absorption of the Dy 
layer on the top. While the precession is equally damped for the Ni thicknesses ranging 
from 7 to 23 nm, an increased damping is found for smaller thicknesses below this. The 
difference in damping of the oscillations is most evident for t^i = 4 and 5 nm. 

The result of the analysis as described before is summarized in figure [5| In this 
graph the data are shown for the samples with the 5nm V, Pd, Dy spin-sink layer and 
the Ni reference. While for the Ni reference, and Ni with adjacent V and Dy layer, 
the frequency dependence is almost equal, indicating similar magnetic properties for 
the different wedge-like shaped samples, the frequency for Pd is found to be somewhat 
higher and starts to drop faster than for the others. The most probable explanation is 
that this difference is due to a slightly different anisotropy for the Ni grown on top of 
Pd in this case. Nevertheless, the magnetic damping found for larger thicknesses ^Ni is 
comparable with the Ni reference. In the upper graph of figure [5] the Gilbert damping 
as a function of the Ni layer thickness is shown. While for the Pd and Dy as a spin 
sink material a additional increase below 10 nm contributing to the damping can be 
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Figure 5. Gilbert damping parameters a and frequency as a function of the nickel 
layer thickness for the intrinsic case and for different damping materials of 5 nm V, 
Pd, and Dy adjacent to the ferromagnet. a is extracted from experiments over a large 
field region. The fits are made using equation 5 and equation 7. In the inset the data 
is shown on a reciprocal scale. Below, the frequency is given (150 mT). The lines are 
guides for the eye. 

For the adjacent V layer, since it is a transition metal with a low spin orbit- 
scattering (light material with low atomic number Z), with a low spin-fiip scattering rate 
and thus a spin diffusion length larger than the thickness ^nm {d <C Agd), no additional 
damping will occur. For Pd and Dy the situation is different: whereas the heavier Pd 
belongs to the transition metals with a strong orbit-scattering (heavy material with 
high atomic number Z), Dy belongs to the rare earth materials. It owns a localized 4f 



Gilbert damping in Nickel thin films 



14 



magnetic moment: therefore, both own a high spin-flip scattering rate and we expect 
the latter two to be in the region where {t ^ Agd). In their cases the thickness of 5nm 
of the spin-sink layer is chosen to be larger than the spin diffusion length (t^M ^ Kd)- 
In this limit the spin current emitted from the magnetic layer through the interface is 



This is called the limit of a perfect spin sink. The additional non-local spin current 
damping is expected to behave inversely proportional with the nickel layer thickness 
^ ^FM- The inset gives the analysis and the data point on a reciprocal scale. The slope 
shows a linear increase for thinner nickel layers, as expected for an inverse proportionality 
for both the Pd and the Dy. Since the value for the intrinsic damping of the nickel fllm 
increases below 4nm this contribution has to be subtracted to reveal the spin-current 
contribution. The value for is then found to be 0.07 for the 2nm Ni/5nmPd fllm, 
which is in the order found by Mizukami by FMR for sputtered Permalloy fllms with 
a Pd spin sink {a^ = 0.04 for 2nmPy/ 5nmPd) [49,50]. A further analysis of the 
thickness dependence of a yields values for the prefactor in [7] for Pd (0.33(3) nm) and Dy 
(0.32(3) nm) with the flt given in the graph. From that value the real part of the interface 
spin mixing conductance in[7|can be calculated. It is found to be = 4.5(5) • 10^^ ficm^ 
for the Ni/Pd and Ni/Dy interface. The increase of the intrinsic damping ai^t has 
been analyzed using an inverse thickness dependence (prefactor of 0.1 nm). While it 
describes the data in the lower thickness range, it can be seen that it does not describe 
the thickness dependence for the thicker range and thus, probably the increase does not 
originate from an interface effect. 

3. 6. Increased damping caused by anisotropy fluctuations: consequences for the 
all- optical approach 

In this last part we want to focus on the deviation from the intrinsic damping a[^t for 
the thin nickel layers itself (^Ni < 4nm). In the low fleld range (10 — 50 mT) small 
magnetization inhomogeneities can build up even when the magnetization appears to 
be still saturated from the hysteresis curve (the saturation flelds are a few mT). For 
these thin layers the magnetization does not align parallel in an externally applied fleld 
any more, but forms ripples. The influence of the ripples on the damping is discussed 
in reference [32]. In the following we adopt this ansatz to the experimental situation 
of the TRMOKE experiment. We deduce a length scale on which the magnetization 
reversal appears for two different Ni thicknesses and relate it to the diameter of our 
probe spot. Lateral magnetic inhomogeneities were studied using Kerr microscopy at 
different applied magnetic flelds [44]. Magnetization reversal takes place at low flelds 
of a -0.5 to 2mT. The resolution of the Kerr microscopy for this thin layer thickness 
does not allow us to see the extent of the ripple effect in the external fleld where the 
increase of a and its strong fleld dependence is observed. However, the domains in the 



totally absorbed within the non-magnetic layer. One can simplify |6] to: 




(8) 
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demagnetized state also mirror local inhomogeneities. For our Nixnm/Si(100) sample 
this is shown in figure [6^) and b) . The domains imaged using Kerr microscopy are 
shown for a 3nm and a 15 nm nickel layer in the demagnetized state. The domains of 
the 15 nm layer are larger than the probe spot diameter of 30 /xm, whereas the domains 
of the 3nm layer are much smaller. 
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Figure 6. a) and b) Kerr microscopy images for the demagnetized state for 15 nm 
and 3 nm. c) and d) corresponding model representing the areas with slightly varying 
anisotropy 



From that observation, the model of local anisotropy fiuctuations known from 
FMR [30,31] is schematically depicted in figure |6]3) and d). A similar idea was also 
given by McMichael [61] and studied using micromagnetic simulations. While for the 
thick film the laser spot probes a region of almost homogeneous magnetization state, for 
the thin layer case the spot averages over many different regions with slightly different 
magnetic properties and their magnetization slightly tilted from the main direction 
averaging over it. The TRMOKE signal determined mirrors an average over the probed 
region. It shows an increased apparent damping a and a smaller resulting from the 
line broadening and different phase in frequency space. While for the thick layer the 
typical scale of the magnetic inhomogeneity is as large as the probe laser spot given 
and only 1-2 regions are averaged, for the thinner film of ^Ni = 3nm many regions 
are averaged within a laser spot, as can be seen in[6|3) and d). Because the magnetic 
inhomogeneity mapping local varying anisotropics becomes more important for smaller 
fields, it also explains the strong field dependence of a observed within that region. 

Figure fTl shows data calculated based on the model, in which the upper curve (i) 
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is calculated from the values extracted from the experimental data for the 10 nm nickel 
layer, curve (ii) is calculated by a superposition of spectra with up to 5% deviation 
from the central frequency at maximum and curve (iii) is calculated by a superposition 
of spectra of 7% deviation from the central frequency at maximum to mimic the line 
broadening. The corresponding amplitudes of the superposed spectra related to different 
Kz values is plotted in the inset of the graph to the given frequencies. The apparent 
damping is increased by 0.01 (for 5%) and reaches the value given in figure ^) for 
the 10 nm film determined for the lowest field values of 30 mT. These effects generally 
become more important for thinner films, since the anisotropy fiuctuations arising from 
thickness variations are larger, as shown by the Kerr images varying on a smaller length 
scale. These fiuctuations can be vice versa determined by the analysis. 




1 ' 1 ' \ ' 1 ' 1 

250 500 750 1000 
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Figure 7. a) Datasets generated by superposing the spectra with the frequency spread 
according to the inset: (i) is calculated from the values extracted from the experimental 
data for the 10 nm Ni layer, (ii) by a superposition of spectra with up to 5% and (iii) is 
calculated by a superposition of spectra owing 7% variation from the central frequency 
at maximum. The average precession amplitude declines faster if a higher spread of 
frequencies (i.e. different anisotropics) are involved. 
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4. Conclusion 

To conclude, we have shown that aU-optical pump-probe experiments are a powerful 
tool to explore magnetization dynamics. Although the optical access to the magnetic 
layer allows an access to the surface only, magnetization dynamics can be explored 
directly in the time domain, resolving different types of spin-wave modes (Kittel mode, 
perpendicular standing spin waves and Damon- Eshbach dipolar surface waves). This is 
in contrast to FMR experiments, where the measured data is a response of the whole 
sample. The obtained data can be similar to the field-pulsed magnetic excitations and 
the Gilbert damping parameter a, needed for the analysis of magnetization dynamics 
and the understanding of microscopic energy dissipation, can be determined from these 
experiments. We have evaluated the contributions of non-local spin current damping 
for V, Pd and Dy. Yet there are limits, as shown for the nickel layer thicknesses 
below 4 nm, where the magnetic damping increases and may overlay other contributions. 
This increase is attributed to an inhomogeneous line broadening arising from a strong 
sensitivity to local anisotropy variations. Further experiments will show whether the 
reduction of the probe spot diameter will improve the results by sensing smaller areas 
and thus reducing the line broadening. 
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